Kondo ground state in a quantum dot with an even number of electrons in a magnetic 

field 
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Kondo conduction has been observed in a quantum dot with an even number of electrons at the 
Triplet-Singlet degeneracy point produced by applying a small magnetic field B orthogonal to the dot 
plane. At a much larger field B ~ £>, , orbital effects induce the reversed transition from the Singlet 
to the Triplet state. We study the newly proposed Kondo behavior at this point. Here the Zeeman 
spin splitting cannot be neglected, what changes the nature of the Kondo coupling. On grounds of 
exact diagonalization results in a dot with cylindrical symmetry, we show that, at odds with what 
happens at the other crossing point, close to _B«, orbital and spin degrees of freedom are "locked 
together", so that the Kondo coupling involves a fictitious spin 1/2 only, which is fully compensated 
by conduction electrons under suitable conditions. In this sense, spin at the dot is fractionalized. We 
derive the scaling equation of the system by means of a nonperturbative variational approach. The 
approach is extended to the B ^ _B*-case and the residual magnetization on the dot is discussed. 

PACS numbers: 71.10.Ay, 72.15.Qm, 73.23.-b, 73.23.Hk, 79.60.Jv, 73.61.-r 



I. INTRODUCTION 

Low temperature transport properties of quantum dots 
(QD) are primarily determined by electron-electron in- 
teractions [0 . When a QD is weakly coupled to the con- 
tacts its charging energy is larger than the thermal en- 
ergy. Therefore, in the absence of a source-drain voltage 
bias, V s d, electrons can be added to the dot in a con- 
trolled way by changing an applied gate voltage, V g H§|. 
Conductance peaks are separated by Coulomb Blockade 
(CB) regions, in which the electron number N is fixed 
(Coulomb oscillations). 

However, there is a well established experimental ev- 
idence for zero-bias anomaly in the conductance in a 
CB region when N is odd This occurs when the 

temperature is further decreased and the strength of the 
coupling to the contacts increased. Such a behavior is 
attributed to the formation of a Kondo-like hybridized 
state between the QD and the contacts |||7j]. Kondo 
effect is not expected when N is even, because of lack 
of spin degeneracy as the QD is in a singlet state. An 
exception is provided when a QD with N even is in a 
triplet state because of Hund's rule. In this case the 
QD shows a very peculiar spin-1 Kondo Effect [p|dTT|. 
The Kondo anomaly is enhanced if the singlet state is 
close to be degenerate with the triplet. Such a regime 
can be achieved experimentally by applying a weak mag- 
netic field orthogonal to the dot, B±, because Hund's 
rule breaks down quite soon. Correspondingly, the QD 
undergoes a triplet-singlet transition, which eventually 
suppresses the Kondo effect. In the following we shall 
refer to the transition above as "TS crossing" . 

If the magnetic field is further increased, orbital ef- 



fects induce the reverse transition from the singlet to the 
triplet state at some "critical field" B± = , in the gen- 
eral landscape of transitions to higher spin states |l4^3[ . 
In this paper we focus on this transition, in the following 
referred to as "ST crossing" . The accidental degeneracy 
between states with different total spin can restore the 
Kondo effect. However, as we will show below, the pres- 
ence of a substantial Zeeman splin splitting (Zss) makes 
the ST crossing very different from the TS crossing. Our 
analysis is an extension of the results of [|14j . 

We consider a vertical QD with contacts located at 
the top and at the bottom of a pillar structure || (see 
Fig. |I]a). Because the confining potential is chosen to be 
parabolic in the radial direction, single particle states on 
the dot are labeled by (n, m, a) where n is an integer, 
m is the orbital angular momentum and a the electron 
spin. We have chosen N = 2 for sake of simplicity, but 
we believe that the pattern we describe is quite frequent. 
For example, it can be extended to the case when A^=4. 

Because of the Zss, the crossing occurring at the ST 
point involves the singlet and the S z = 1-component of 
the triplet state only. The analysis of the quantum num- 
bers of the N ± 1 states visited by virtual transitions (co- 
tunneling processes ) shows that orbital and spin quan- 
tum numbers of the electrons involved are "locked to- 
gether" (Fig. 2a). Tunneling from the contacts does not 
conserve n, but, since the geometry is cylindrical, it con- 
serves m and the spin. In particular, if the dot is in the 
triplet state, only (m = 0, J,) -electrons can enter it, while 
is it is in the singlet state, then only (m = 1, t)-electrons 
enter. This envisages a one-channel spin 1/2-like Kondo 
coupling different from what occurs at the TS point. 

The existence of a Kondo effect between states belong- 
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ing to different representations of the total spin is very 
peculiar. Being the total dot spin of both the degener- 
ate levels integer (S — 0,1), spin compensation has to 
be incomplete. In this work we show that it is as if the 
total spin were decomposed into two fictitious spins \ 
which we refer to as S r and S e g p5| . We stress that 
these dynamical variables involve both orbital and spin 
degrees of freedom which are locked together. The mag- 
netic field B* favors antiferromagnetic (AF) coupling be- 
tween the delocalizcd electrons of the contacts and S e jj, 
which quenches it in the strong coupling limit. S r is the 
residual spin i at the dot 0. CB pins N to an even 
value, but the correlated state that sets in will have half- 
odd spin, i.e., "fractionalization of the spin" at the QD 
site will take place [0. B* only affects the magnetiza- 
tion related to S r , but we will not discuss such a point in 
this paper. On the contrary, detuning B from B* corre- 
sponds to an effective magnetic field SB — B — B*, which 
affects Kondo screening. 

On the other hand, at temperatures T > Tk ( Tk 
being the Kondo temperature), there is no Kondo corre- 
lation. Furthermore, being the dot in CB, the screening 
by the delocalized electrons is quite negligible. As long 
as B < B*, the total spin is 0. For B > the total spin 
is 1. The spin density on the dot has a drastic change 
at while the charge density is only slightly affected. 
This is shown in Fig. 3a: for B < B„ (singlet state) the 
spin density is zero, while it is half the charge density for 
B > B* (triplet state S = 1 with z— component S z = 1). 

This shows that the total spin of the dot is the relevant 
dynamical variable in the correlation process that sets in 
when temperature is lowered, while the charge degree of 
freedom is frozen for B ~ £>*. 

In Section IV we set up a variational technique similar 
in spirit to Yosida's approach |M for the construction of 
a trial wavefunction of the correlated Ground State (GS). 
We start from a Fermi sea (FS) of delocalized electrons of 
the contacts, times the impurity spin wavefunction and 
we project out states which do not form a singlet at the 
dot site. This is a modified version of the Gutzwillcr 
projector and the results we obtain qualitatively repro- 
duce the strongly-coupled regime of the Kondo model 
[Rl9H2lj. Our formalism also allows us to discuss the 
SB 0-case and to calculate "off-critical" quantities, like 
the magnetization, as a function of SB. In H] we esti- 
mated the Kondo temperature Tk- It is slightly below 
the border of what is reachable nowadays within a trans- 
port measurement (Tk ~ tens of mK). Instead, an Elec- 
tron Paramagnetic Resonance (EPR) experiment could 
provide evidence for the fractionalization of the spin. At 
T < Tk, the system goes across the Kondo state when 
B is moved across B* . Partial compensation of the dot 
spin due to the local screening implied by the Abrikosov- 
Suhl resonance takes place at the QD site. The residual 
magnetization is the one of a spin- 1/2 in all the range 
B ~ B*. Therefore, in lowering the temperature, the 



response to the microwave field should also appear all 
the way from B < £>* to B > B* and a temperature- 
dependent Knight shift of the frequency should occur. 
Moreover, this feature should be detectable, due to the 
enhanced susceptibility of the Kondo state. 

Differently from the scenario described here, the TS 
crossing has a nonsymmetric behavior across the degen- 
eracy point. Spin 1 Kondo effect should take place for B 
prior to the degeneracy point, It is renforced at the TS 
transition, but no Kondo effect at all should occur when 
the GS is a singlet. 

The paper is organized as follows: 

• In Section II we show that B* makes only four levels 
involved in the Kondo coupling, in analogy to the 
single impurity spin i Anderson model (AM). We 
discuss the strongly coupled regime of the system 
by briefly recalling the scaling theory of the AM and 
specify what we mean as spin fractionalization. 

• In Section III we employ the Schrieffer- Wolff trans- 
formation to map our system onto an Effective 
Kondo model. We show that the correspond- 
ing interaction between S e ff and the spin of the 
electrons from the Fermi sea is antiferromagnetic 
(AF), which generates the strongly coupled Kondo 
regime. 

• In Section IV we construct a trial GS using a mod- 
ified version of the Gutzwillcr projection method. 
We calculate the energy, which depends on the ex- 
change interaction: j. We show that j scales ac- 



cording to Anderson's poor man's scaling 22 



• In Section V we extend our technique to the SB ^ 
0-case. Because of the presence of SB, a second 
scaling parameter, h, arises. We compute the en- 
ergy and the magnetization as a function of h and 
discuss how this is related to the bare field SB. 

• In Section VI we report our final remarks and con- 
clusions. 



II. THE MODEL FOR THE DOT AT THE LEVEL 
CROSSING (B = B,). 

Using exact diagonalization , two of us studied few 
interacting electrons in two dimensions, confined by a 
parabolic potential in the presence of a magnetic field B 
orthogonal to the dot disk (z— direction). We chose the 
energy separation between the levels of the confining po- 
tential hwd to be of the same size as the effective strength 
of the Coulomb repulsion U (~ 4meV), which makes the 
correlations among the electrons relevant. The magnetic 
field favors the increase of the total angular momentum 
in the z~ direction, M, as well as of the total spin of the 
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dot, S, in order to minimize the Coulomb repulsion. The 
independence of the spectrum of the QD energy levels 
can be monitored via a linear transport measurement, 
by attaching electrodes to it in a pillar structure (Fig. [I]) 
||,|]J|]. At special values B = B*, one can identify sev- 
eral level crossings between states with different M and 
S. 

In this paper we shall focus onto a QD with N = 2 
electrons, but the features we discuss here are appropri- 
ate for dots with any even N. The parameters of the 
device can be chosen such that B* is rather large, which 
lifts the spin degeneracy, because of Zeeman spin split- 
ting (23) . We define -B* as the field at which the GS of the 
isolated dot switches from singlet ( 2 Sq; S = 0, M = 0) 
to the triplet state lowest in energy (the state S z = 1 of 
the three states 2 T^ ±1 : S=1,M = 1 and S z = 0,±1). 

Only two N — 2 states are primarily involved at the 
level crossing, namely, the 5 = singlet and the S = 
1, M = 1,S Z = 1 component of the triplet (see Fig. |2^l). 
At the same B, the low lying N = 1 and N = 3-GS's 

are spin doublets. We shall denote them as 1 D\ and 

2 

D\, respectively. Lifting the spin degeneracy makes the 

2 

S z = I components lower in energy in both cases. If V g 
is tuned in the CB region between N = 2 and N = 3 
(Fig. |l|b), only the four states listed above are primarily 
relevant for the conduction. 

We now construct a model by taking the 1-particle 

state , 1 D°i, as the "vacuum" of the truncated Hilbert 

2 

space S (which we denote by |0)). Let us take |0) to be 
the initial state of the QD. The dot can exchange one 
electron at a time with the contacts. A (q, m, f )-electron 
entering the QD generates a transition between |0) and 
2 T 1 1 . Alternatively, the lowest energy state available for a 
(q, m, |)-electron entering the QD is the singlet 2 5q (here 
q is the modulus of the momentum of the particles and 
m is the azimuthal quantum number, which is conserved 
by the interaction with the dot) . It is useful to define the 
operators d\ 2 which create the two N = 2 many body 
states I a) (a = 1, 2), belonging to S, by acting on |0): 

|1)^| 2 S °)=4|0); |2> = | 2 I?>=4|0) (1) 

The operator d\/d\ adds a (m = 0, |)/(m = 1,T) elec- 
tron to the dot. Finally, the state |3) is defined as: 
|3> = I 3 D\) = dld\\0). The operators d 12 defined 
above cannot be associated to any single particle func- 
tion on the dot. Nevertheless, as far as low-energy exci- 
tations only are involved, it is a good approximation to 
take them to be anticommuting. To show this we refer to 
Fig. |], where we report the calculated quasiparticle spec- 
tral weight, when an extra electron is added to the dot 
in one of the two possible GS's at B = B* and TV = 2. In 
Fig. |^a we show the weight for addition of an extra par- 
ticle to the state | 2 Tl) to give a state with total angular 
momentum M = 1. This is close to 1 for a |-spin but 



is practically zero at low energies for an f spin because 
states with S — 3/2 are higher in energies. According 
to the definitions of eq. ([!]), this implies that, while the 
state djdJlO) can be normalized, one has (<4) 2 |0) ~ 0. 
Analogously, Fig. ^b shows that there is an energy shift 
of the spectral weight peak when a f or a J, spin particle 
is added to the state | 2 5q). This shift corresponds to 

the Zeeman spin splitting between the states 3 D\ and 

2 

3 D l _i. We ignore the latter that is placed at higher en- 

2 

ergy. In Section III we argue that inclusion of this state 
does not qualitatively change our results. Then, within 
this approximation that includes the lowest energy state 
only, we take (d{) 2 |0) w 0. So, we are alleged to assume 
(d\ 2 ) 2 = 0. Hence, the d-operators behave as they were 
single-particle fermion operators, although, in fact, they 
create many-particle states. 

We write the Hamiltonian for the QD as: 

Hd = (-d(n\ + n 2 ) + Unin 2 (2) 

where n a = d* a d a , £d is the energy of the degenerate 
N = 2 levels and U is the charging energy for adding a 
third electron. 

We assume the electrons in the leads to be noninter- 
acting. Let bj qa be the annihilation operators for one 
electron in the left (j — L) and the right (J = R) con- 
tact, respectively, and let e q be the corresponding energy 
(independent of j and a for simplicity) . The Hamiltonian 
for the leads is: 

qc j=L,R 

We denote by 7i = Hjj + Hi the sum of the two Hamil- 
tonian terms given by eqs.(^,||). They do not change the 
particle number N on the dot. Instead, N is changed 
by tunneling between the leads and the dot, according to 
the above mentioned selection rules. Let and be 
the amplitudes for single-electron tunneling between the 
dot and the R or L contact, respectively. Removal of one 
electron from the contacts is described by the following 
linear combination of the 6j gcr 's: 

= /in 12 , it. If ( Fl hhqa + TRbR 4° ) ( 4 ) 
V\ T l\ + \ T R\ 

For a point-like QD, there will be no tunneling term 
involving the combination of bjka operators linearly in- 
dependent of 4> qa - The terms to be added are (a — 1,2): 

V a0 = <o(<?)4<V Q (5) 

for N = 1 -► N = 2 and: 

V?,a = Y V3a O^a'Vc, j (6) 

Q 
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for N = 2 — > N = 3(a ^ a) . Because spin is conserved 
in the tunneling is <j\ and 02 =T an d a ^ a. Her- 
mitian conjugate operators of V a o and V^a decrease the 
electron number on the dot. 

This model Hamiltonian has the form of the nondcgcn- 
erate Anderson Model (AM). In |0| we have considered 
a special choice of the V g fine tuning (see Fig. |l|b) , which 
corresponds to the symmetric case: the singly occupied 
localized impurity level has energy e<j = — £7/2, while the 
doubly occupied one has energy td + U = U/2. Also, 
Tj (j = L, R) have been taken independent on j. 

Integration over the contact fields 4> qa in the partition 
function leads to the constrained dynamics of the dot be- 
tween the two degenerate states \a), labeled by a = 1, 2 
and coupled to an external fluctuating field X(r) @. 
As shown in fl4f| , the main contributions to the partition 
function are the quantum fluctuations of an effective fic- 
titious spin i S e ff, due to the coupling among the dot 
and the leads. This spin flips from 1 to 2 because of 
cotunneling processes in which an (to = 1, j) electron is 
added/removed in the dot while an (to = 0, J.) electron 
is removed/added (see Fig. 2a ) Its quantum dynam- 
ics produces a correlated state between dot and leads in 
which the charge excitations, being higher in energy, are 
fully decoupled and suppressed. Four states are available 
prior to correlation, corresponding to S = 0,1, two of 
which are degenerate. They factorize according to the 
decomposition law: 



S — Sf ~\~ s, 



eJJ- 



(7) 



where S r is a residual spin 1/2 whose wavefunction fac- 
torizes close to -B* . Correlation with the contacts is gov- 
erned by a single channel spin ^ AF Kondo Hamiltonian 
for S e ff, which is derived in Section III and discussed in 
detail in Appendix A. In the corresponding GS, S e g is 
fully compensated as in the usual Kondo Effect at zero 
magnetic field. 

In Sect. Ill and IV we show that, provided that the 
coupling to the leads is symmetrical for a—\,2, the GS 
is doubly degenerate at B = (the limiting states are 
derived from the \g x ) (A = ±) of eq.(|l6|), having en- 
ergy e(£,j) given by eq.(ptj|)). The degeneracy disap- 
pears when is SB = B — B* ^ (see Fig. 2b). Such 
a spin fractionalization closely resembles quantum num- 
ber fractionalization in strongly-correlated 1-dimcnsional 
electron systems 129] . 



III. EFFECTIVE AF KONDO INTERACTION. 

In the previous Section we have introduced the four 
states \i), i — 0, 1,2,3 , which are primarily relevant to 
the low-temperature dynamics of the QD (see eq.(Q)) 



Tuning V g as described in Section II (see Fig. |l|b) allows 
only the N = 2 states, |1) and |2), to be the initial and 



final states of any process. Let S be the subspace they 
span. The N ^ 2-states, |0) and |3), play a role only 
as virtual states allowing for higher-order tunneling pro- 
cesses. Then, it is desirable to describe the dynamics of 
the dot in terms of an effective Hamiltonian acting on S 
only. This is accomplished in this Section, by employing 
the Schrieffer- Wolff (SW) transformation at second order 
in the tunneling amplitudes. 

According to eqs.(||,|| ) the elements of Tio are diagonal 
in the Fock space and are given by: 

Hoo = Hi; Hn — H22 = H[ + e^; H 3i = H[ + U + 2ed 

(8) 

Hq3 — H12 — 0, while the tunneling part of the Hamil- 
tonian of eq.s ((^^) provides the remaining off-diagonal 
terms i?y (i, j — 0, 1, 2, 3). 

The effective Schrodinger equation projected onto S is: 

P[(H - E) - yt(i _ p)(H - S)- X F]P|*) = (9) 

where P is the corresponding projection operator. 

For the sake of simplicity we shall take v^, a (q) and 
v a o{q) to be independent of a. Then, to second order 
in v, eq.(|| ) becomes (see Appendix A for details): 



qq' 



qq' a 



(10) 



where a^ x , are the spin ^ matrices (a = x, y, z) and 
Seff is a spin-i representation of the rotation within S, 
which is well defined, provided the constraint m +ri2 = 1 
is satisfied. Indeed, in this case, (S e ff) 2 = f (ni + ri2 — 
2n\ri2) — |. More explicitly: S z e g— ^{d\d,2 — d\di) and 

the ladder operator is S~g = d\d2- 
The interaction matrix elements are 



Q q 



C-,1 



U-6 



T , — 2 



^3(9)^3(9') _ vo(q)vo(q 

e d + U - e„ 



(11) 



corresponding to the potential scattering term and the 
Kondo coupling term, respectively. 

The potential scattering term Q qq i provides a one-body 
interaction in the "charge" channel. It is irrelevant for 
the physics of the Kondo effect, which is related to inter- 
actions in the "spin" channel [^Tj| , and we shall disregard 
it in the following. 

By taking q at the Fermi surface, the Kondo term in 
eq. (|Ty) becomes: 



H K = H cS - Hi 



q'X'- 



(12) 



</</ 
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Tuning V g such that the chemical potential of the leads 
is the zero of the single-particle energies, the Kondo cou- 
pling strength in eq.fO) takes the usual form: 



J=\Vf 



1 



1 



(d + U e d 



(13) 



ed Ri —U/2 implies J > 0, i.e., the interaction between 
the spin density of the Fermi sea at the location of the 
dot and the effective spin S e ff, described by eq. (p^) , is 
antiferromagnetic. Thus, the low-energy physics of the 
system will be controlled by the strongly coupled AF- 
S = i-Kondo fixed point, where a spin singlet takes place 
at the impurity site. Because of the way S is related to 
SeS (eq.(0 )), total compensation of S e g leads to a spin 
S = 5 at the QD. 

In the next Section we will make use of such an equiv- 
alence to construct a formalism for the strongly coupled 
fixed point of our system. 



U, where higher-energy states are "projected out" by the 
interaction and the GP method works successfully p7| . 
Our approach is similar in spirit to Yosida's variational 
technique |Tg| ]. However, here we are mostly interested 
in the "macroscopic" variable <r(0). Hence, it is the only 
lead operator involved in the construction of the vari- 
ational state. This makes our technique much simpler 
than the one in jl8| , because our approach requires just 
one variational parameter for the trial state. 

Here we list the basic assumptions concerning the 
model Hamiltonian and the trial variational state. We in- 
troduce our approach for the simple spin-i Kondo Hamil- 
tonian (which corresponds to B = £>»). In Section V we 
generalize it to the SB ^ 0-case. To simplify the nota- 
tion, we drop the suffix e jj from S e g all throughout this 
Section and the next one. 

a. The Hamiltonian is (see eq.(fl2|)): 

H = H i + H K = H i + Jaa(0) ■ S (14) 



IV. GUTZWILLER PROJECTION. 

In Section II we have shown that cotunneling at the 
ST point strongly involves just two among the four avail- 
able states of the isolated dot, which can be described 
in terms of S e jj, including both orbital and spin degrees 
of freedom locked together. Instead, the other relevant 
degrees of freedom of the dot can be lumped into another 
pseudospin variable S r (S r = |) which decouples below 
T K , because the dot is tuned at CB (eq.(^)). According 
to Section III the dynamics of correlations involves the 
true spin density of the delocalized electrons at the dot 
site ct(0). Therefore, the study of the magnetization in 
the correlated state requires the knowledge of <?(0). 

In Section III we discuss the assumptions under which 
the isotropic single channel spin | AF Kondo model of 
eq.(|lj) describes the physics at the fixed point. The cor- 
related state is a Nozieres local Fermi Liquid (FL) [p6| , 
with a spin singlet at the origin. In our case there is a 
substantial magnetic field B ~ £>» and a small detuning 
SB ^ is likely to occur. 

In this Section, we construct an approximated GS of 
the strongly coupled system, described by H e ff (eq.(^2[)). 
In order to describe the scaling of S e g toward total com- 
pensation, we use a variational method based on the 
Gutzwiller projection (GP) technique. This qualitatively 
reproduces the main features of the correlated singlet 
state. Detuning B off £?* allows for probing the QD 
magnetization and for studying the magnetic response 
in the correlated state. So, in Section V we extend our 
approach to the B ^ B^-c&se. 

As J increases, states other than a singlet at the im- 
purity become higher in energy and can be "projected 
out" from the physical Hilbert space. This is quite sim- 
ilar to what happens in the 1-d Hubbard model at large 



The impurity is located at x — 0. The scattering 
is s-wave, so that the model is effectively 1-d and 
can be defined on a lattice with spacing a and pe- 
riodic boundary conditions at L — N\a. (N\ is the 
number of lattice sites ) |2j| . Hk provides an effec- 
tive interaction among the delocalized electrons. At 
each stage of the scaling process, we take the Slater 
determinant \FS) = JX II 



q<0 u qa 



|0) as a reference 

state (FS: Fermi sea). The c qa operators annihi- 
late the quasiparticles and |0) is the quasiparticle 
vacuum. The spectrum is linearized around the 
Fermi point qf- D is the bandwidth. Therefore, 
the lead Hamiltonian is Hi = J2 q a e q<? c \<7 c qv j where 
£(?ct = vfI- Here vf is the Fermi velocity and q is 
the momentum measured with respect to the Fermi 
momentum. The bandwidth is 2D = 2-kvf/ol- The 
density of states at the Fermi level is ^(0) = Ni/D 
and is assumed to be constant during the scaling. 

b. One can represent the real space field operator in 
terms of the c-operators: 



(15) 



It satisfies the anticommutation relations 
{fc^( S ')} = i* aiB /^. 

The quantity an a (0) = a^)J(0)^> ff (0) plays the 
role of the occupation number at the impurity 
site. This is consistent with the operator identity 
(an CT (0)) 2 = an (T (0). The corresponding spin am- 
plitude due to the delocalized electrons at x = is 
a<?(0) = ^ai/jt{Q)v aa 'i' a ' (0). The average value of 
(a<r(0)) 2 on \FS) is: (FS\(aa(0)) 2 \FS)=3/8. This 
is a consequence of averaging over configurations 
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with zero or double occupancy at x — (total spin 
0) and configurations with single occupancy at the 
same point (total spin 1/2). 

c. The correlation between the impurity and the FL is 
accounted for by projecting out of the uncorrelated 
state the components with a triplet or a doublet of 
the total spin at the impurity site, S + aa(0). 

Denoting the projector as P g , the variational trial 
state is defined as: 



\g x )=P g \FS,\)^P g (\FS)®\\)) 



(16) 



Here |A) is an eigenstate of S z with eigenvalue j\ = 
±i and g is the variational parameter. The explicit 
form of the Gutzwiller projection operator P g will 
be given below. So far, it is enough to say that all 
the components of the state other than a localized 
singlet at x = are projected out, as g changes 
from zero to the limiting value, g = 4/3. 

d. At B = B*, we assume the usual one parameter scal- 

ing process when D — ► D — SD. g approaches 
4/3, and, in order to to guarantee that at each 
step the total energy is at a minimum, the pa- 
rameter v(0) J is correspondingly renormalized. If 
SB = B — £?* ^ 0, a second parameter appears in 
the model, related to SB. We shall refer to such a 
parameter as h and will introduce it in the next Sec- 
tion. It is generated by the shift in the quasiparticle 
energies, according to their spin and will determine 
the "off-critical" magnetization (see point e). 

e. According to eq.(Q), the effective magnetization of the 

dot is given by the average of the total spin at x = 
on the trial state \g x ) (shortly denoted by 



}g,8B) 



lim {M d (g, SB) - M d {g, SB = 0)} 

3^4/3 



M, 



aa z (0)) g:SB ] 
imp + m(h) 



(17) 



Here we have temporarily restored the suffix e ff, 
for clarity. The dot magnetization M d includes the 
local moment on the impurity site Mi mp , and the 
induced magnetic moment of the delocalized elec- 
trons from the contacts, m(h). 

At SB = both Mi mp and m(h) vanish, due to the 
compensation of the effective spin S e ff of the stan- 
dard spin | AF Kondo problem. Md{g, SB = 0) is 
the residual magnetization at B = B* associated 
to the expectation value of S r of eq. (0) and is not 
discussed in this paper. However, the spin suscep- 
tibility of S r is expected to be much smaller than 
the Kondo susceptibility of S e g. 



One remark concerns the fact that the giromag- 
netic factor in eq.(p7|) might not be the bare one, 
appearing in the electron Zeeman spin splitting. 

In Section V we work out the form of Mi mp and 
m(h) in the strongly coupled scaling regime. 

Let us now construct P g . 

According to Noziere's FL picture of the correlated 
Kondo state [^6), the main effect of the interaction be- 
tween the FS and the impurity taken at the AF fixed 
point is to fully project out the components of the trial 
state other than a singlet of the total spin at the impurity 
site. Let us consider the operator 



(aa(0)f 



(18) 



It will give zero when acting on a state that has a doublet 
on the impurity site, (i.e. ((acr(0)) 2 ) = 0), —3/4 when 
acting on a singlet, 1/4 when acting on a triplet. Among 
the three possibilities listed above, the state correspond- 
ing to the triplet is the one highest in energy. We fully 
exclude it from the trial state from the outset and grad- 
ually project out the doublet by varying the parameter g 
from to 4/3. Accordingly, we define P g as: 

P g = (l--g- gaa{0) ■§){!- 40^(0) • S) = 



1 - -g + g (aa(0)) 2 ) -W(0)-5 (19) 



This operator projects out the components other than 
the localized singlet. When g=4/3 the projection is com- 
plete. 

In this section we report the result for the zero mag- 
netic field case. The key points of the calculations are 
summarized in Appendix B. 

The expectation value of the Hamiltonian in eq.jl^) 
differs from the energy of the reference Fermi sea \FS), 
E n = —2nvpNi/a, just by a term of order O(l) in the 
particle number iVj. In units of D, this correction is: 



_ i / VW) 



2 i±il^ (20) 



where N[£] is norm of the trial state (see eq.(Bl)) and 
we have introduced the dimensionless variables j = 
3i>>(0) J/(8Ni) and £ = |(| — g) for convenience. 

It is worth stressing that e(£,j) does not depend on 
the polarization of the impurity, A. Such a degeneracy 
disappears when SB ^ |29|]. 

For each j oc J/D, eq.(|20| ) takes a minimum as a 
function of £ (i.e. g). The minimum £(j) w.r.t. j is given 
by: 



f + 1 



(21) 
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Both the triplet and the doublet are fully projected out 
of the trial GS when g — > 4/3, i.e. when £ flows to zero. 
Eq.(pTl) shows that this limit corresponds to j — > oo. 

The derivative of eq.( |2l|) w.r.t. £, reproduces the poor 
man's scaling law for j ]22{ provided g is small (that is, 
£ is close to 1/2) : 



d(ln£) 



(22) 



Eq.(p2[) has been worked out in |2^] within perturbation 
theory, provided g is small, that is, £ is close to \. In our 

!) holds all the 



simplified approach , eq.(21) (not eq.(|2 
way down to the fixed point. 



V. THE MAGNETIC MOMENT 

In this Section we qualitatively discuss what happens 
to the energy and to the magnetic moment of the dot 
when a magnetic field is present and B is possibly de- 
tuned off _B*. 

Here we do not consider the in detail the magnetic mo- 
ment associated to S r . Instead, we focus on the one in 
the third line of eq.(|l7|). We generalize the construction 
of Sect IV to the SB ^ 0-case by shifting the energies 
of the quasiparticles of the FS according to their spin. 
Two corrections to the magnetic moment arise due to 
this shift: a) {S e g)g^B may not scale down to zero, 6) an 
induced moment of the delocalized electrons m{h) may 
arise. We first discuss the latter contribution. In the ab- 
sence of correlations, because B„ is the degeneracy point, 
the energy difference of electrons involved in cotunneling 
at the Fermi level is the same for both spins. There- 
fore, provided also the density of states of the conduction 
electrons is the same at the Fermi energy for both spins, 
B* will only be responsible for changes in the tunnel- 
ing Hamiltonian due to orbital effects [pOt , which do not 
influence the flow to strong coupling in a relevant way. 

This is not the case if SB ^ 0. Qualitatively, there will 
be changes in the number of the conduction electrons 
involved, according to their spin a = ±^ as follows: 



Ni/2 



Ni 



SB 

r 

D 



(23) 



The interaction modifies eq.(|23|) by scaling, so, we gen- 
eralize it by defining a dimcnsionlcss quantity h in terms 
of the expectation value of the quasiparticle spin at the 
origin evaluated on the reference state \FS): 



(FS,X\aa z (0)\FS, X) 5B ^o = 
(FS,X\\ (an T (0) - an;(0)) \FS, A) = \h x 



(24) 



(h\ = A ■ h, where A = ± is the polarization of the impu- 
rity ). Eq.(|24]) states that the magnetic field produces a 



spin density at the impurity site proportional to h. When 
SB = 0, rotational symmetry of the FS implies h = un- 
less a local paramagnetism takes place, so that h is a 
homogeneous function of SB. 

Eq. ( p4] ) is enough to generalize eq. ( |2Q ) to the SB ^= 0- 
case. Using the results of of Appendix B, we get: 



[2(1 - g + e) + Ml - 2Q](1 - h x ) - 2j(l + h x ) 
[(l + hx) + ?{l-h x )] 



(25) 



As in the previous Section, the minimum of cq.(|25|) 

w.r.t. £ (de/d£, = 0) fixes - (1 + h)/2j. The 

states |g A )(A = ±) are no longer degenerate, except for 
h = and they represent the ground and excited state, 
depending on the sign of h ( Fig. 2b). At h = the 
delocalized electrons do not have an intrinsic magnetic 
moment close to the fixed point: m(h = 0) = 0. Because 
of the degeneracy, the magnetization at h ~ is obtained 
from the average E/J = (e + + e_)/2. Minimization with 
respect to £ yields the leading large j-correction to the 
energy for h ^ 0, due to the delocalized electrons of the 
leads: 



E 
1 



1 

jo 



jo 



0(h 4 



(26) 



Hence, the contribution from the delocalized electrons to 
the magnetic moment is linear in h as expected from the 
exact result ^cj . However, the relation between h and SB 
cannot be established within this model, unless we state 
some direct link between the strong coupling fixed point 
and the uncorrelated state previous to scaling. This link 
can be inferred from the educated guess that the response 
of the conduction electrons to SB, given by eq.(^) in the 
absence of correlations, smoothly evolves into the one of 
the correlated state, m(h), in the form: 



m{h) = lim (g x \aa z {0)\g x ) j!h 



(here Tk is the Kondo temperature). 

Indeed, by using the insight we have from the exact 
solution, we infer that, the fraction of delocalized elec- 
trons involved in the correlated state is expected to be 
0(Tk/J)- Because Tk is a scale invariant, both h and j 
scale to larger values in the same way as £ — * 0. 

Then, by comparing eq.(p7[), supplemented with the 
derivative of eq. (p6[) , to the magnetization of the uncor- 
related state, |(iV T - Ni)/Ni = g*n B SB/4D, we obtain: 



1 dE 
2J~dh 



TTF^T ~ MBS 



SB_ 



(28) 



Thus, we recover the expected result that the spin sus- 
ceptibility is cx 1/T K 111- 

This concludes our discussion of m(h). 
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According to eq.([L7|), the total magnetization also in- 
cludes the term Mi mp . This term just provides a sub- 
leading correction to the total magnetization. 

Indeed, Mi mp is given by 

1 £ 2 
M lmp (t h\) = -^ B g*\ TT ^—— , (29) 

l-hx + ? 

i.e. this correction to the fractional spin of the dot de- 
creases as the second power of £ when £ decreases in the 
strong coupling limit. 

VI. CONCLUDING REMARKS 

Kondo coupling is the striking realization of non- 
perturbative tunneling between a QD in the CB region 
and contacts. It requires the degeneracy of the GS dot 
level, what certainly takes place when N is odd. In the 
even-iV case, the GS can be a triplet, because of Hund's 
rule and corresponding partial screening of the 5 = 1 
spin on the dot has been observed Ji],^]. The proximity 
to a Triplet-Singlet (T-S) transition enhances the cou- 
pling and increases the Kondo temperature M| . A small 
magnetic field driving the system across the level crossing 
can tune the intensity of the effect due to the T-S con- 
version. This enhanced Kondo conductance was indeed 
measured for N = 6. 

In a different mechanism for the Kondo effect 

in an even- A QD has been proposed, considering states 
belonging to different representations of the total spin, 
5=0,1, at the accidental crossing between the singlet and 
the triplet state lowest in energy. In |l7| B is taken in 
the plane of the dot and the crossing is attributed to the 
Zeeman spin splitting term. Because the latter is usually 
very small, it is unlikely that this can occur anywhere else 
but at the point mentioned above. Then the mechanism 
discussed by |||l0| seems to be more likely |H] . 

On the contrary, by increasing Bj_ up to a suitable 
value B* , the degeneracy of the triplet state is lifted and 
the singlet state ceases to be the GS of the system. In JlJ] 
we estimate quantitatively the parameter values of the 
optimal device by using results of exact diagonalization 
for few interacting electrons in a dot |23| . A vertical QD 
with a magnetic field along the z— axis, undergoes a S-T 
transition mainly due to the Zeeman term, which favors 
larger angular momenta along z, and to the correlation. 
The spins of the electrons in the contacts are still unpo- 
larized at the Fermi energy if the contacts are bulk (i.e. 
3-dimcnsional) metals. This allows Kondo coupling to oc- 
cur at the crossing between the two lowest levels among 
the four with S = 0, 1. In |l7j the zero voltage anomaly in 
the differential conductance is discussed by including just 
these two levels. The fact that two different spin states 
are involved is reflected in the asymmetry in the trans- 
mission probabilities w.r.t. the spin index. This can be 



traced back to a term which couples the charge and spin 
degrees of freedom in the exchange Hamiltonian that can 
be derived from a Schrieffer- Wolff transformation. How- 
ever, this term is irrelevant as the system flows toward 
the strongly-coupled fixed point, consistently with the 
full decoupling of the charge and the spin degrees of free- 
dom in such a limit, due to the widely different energy 
scales for the two excitations |2^|. The Schrieffer- Wolff 
transformation (see Section III) maps the problem onto 
the AF spin \— Kondo Hamiltonian. Coupling is better 
described by splitting the total dot spin S = 0, 1 into two 
spins |, which we refer to as S r and S e g (see eq.(^)). S e g 
is totally screened, independently on while a residual 
spin | survives, whose magnetic moment depends on B* . 

In this work we have concentrated on the magneti- 
zation at the strong coupling fixed point. A variational 
approach is used for the GS, based on the Gutzwiller pro- 
jection method, which, as far as we know, has never been 
applied to the Kondo problem. This technique allows for 
studying a small detuning of B off £?* and for discussing 
the dependence of the magnetic moment on SB = B—B*. 
The basic assumption is that SB ^ does not move the 
system away from the strongly coupled state. The scale 
of magnetic field at which the Kondo state is disrupted, is 
expected to be quite large (it was estimated in to be 
of the order of 1 Tesla). A similar conjecture in the case 
of magnetic impurities in metals has been formulated by 
Y. Ovchinnikov et al. ]32]], starting from a perturbative 
approach in a reduced Hilbert space. 

When SB ^ 0, besides the Kondo coupling parameter 
j oc v(0)J, our variational correlated state includes a 
new coupling parameter which is related to the detuning 
magnetic field, h\ oc SB/Tk- As the system flows toward 
strong coupling, this parameter increases together with 
j. Within our formalism we show that j scales according 
to Anderson's poor man's scaling, provided the ration 
h/j keeps small (< 1). 

We have proved that the magnetic moment on the dot 
has a term linear in h (see eq.s( |l7| , ^7|) ), which arises 
from the screening due to the electron density of the 
Abrikosov-Suhl resonance located at the dot site. Still, 
more work has to be done in order to provide the relation 
between h and SB (in our formalism h is introduced as 
the shift of the quasiparticle energies corresponding to 
u{Q)SB). 

The QD GS includes what we referred to as a spinon, 
carrying a magnetic moment S r = 1/2 but no charge. 
The splitting in energy with magnetic field arising from 
this term will depend on the spinon wavefunction itself. 
At this stage, an estimate of the renormalized giromag- 
netic factor is impossible. However, the residual half in- 
teger spin on the dot, together with the screening effects 
in the GS at non zero SB, has important experimental 
implications. 

In |l4| Tk has been estimated to be rather low, i.e., 
of the order of tens of mK. At such a low temperature a 
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transport measurement is cumbersome. We proposed to 
probe the fractional magnetization of the QD by means of 
a magnetic resonance experiment with microwaves. The 
energy splitting of the dot magnetic moment can be de- 
tected by observing resonant absorption as B is slightly 
detuned from B*. 

A Knight shift should be measured, depending on the 
sign of SB, which is a consequence of the term m(h) 
in eq.(^). In Fig. || we plot the radial behavior of the 
z— component of the total spin in the GS when B is 
slightly less or larger than B*. The dot is isolated and 
in the Coulomb blockade regime at N = 2. In absence of 
Kondo coupling, screening from the delocalized electrons 
of the contacts does not occur and the magnetic mo- 
ment changes abruptly from zero to its maximum value 
at the transition. On the contrary, in the Kondo regime, 
the spin susceptibility of the conduction electrons on the 
dot is expected to be large and a temperature dependent 
screening should be present in the neighborhood of 
This gives rise to a smooth crossover in the paramagnetic 
resonance close to -B* at low temperature. 

In conclusion, we have discussed the new possibility 
that Kondo effect arises in a QD with N = 2 at low 
temperature, when the dot is tuned to the ST point and 
the coupling to the leads increases. This happens when 
the GS of the unperturbed dot changes from 2 Sq (B < 
£?*) to 2 T\ (B > B*), due to level crossing. 

The theory of |l4| and of this work is not applicable 
to the TS crossing. Orbital effects are dominant in our 
case in producing the crossing as well as the properties 
of the many body states with N + 1 or N — 1 particles, 
as described in the text. All our results about the nature 
of the states involved are totally independent on the Zss. 
However, at the ST point the Zeeman spin splitting is 
quite substantial, and favors a certain component of S z . 
Its role becomes important in two respects: 

a) Only two A-electron levels cross at the TS point 
(not four oh them, as is the case of the TS point), which 
we label by 1 and 2. 

b) states at N, N + 1 and N — 1, although they 
are many-body levels, have quantum numbers such that 
there is just one channel by which an electron can be 
virtually added or subtracted to the dot in the A-state, 
because spin and orbital momentum of the extra electron 
are strictly "locked together" in the selection rules. 

The spin of the QD becomes \ without changing the 
average occupancy We propose to probe spin frac- 
tionalization with a magnetic resonance experiment. A 
continuous shift in the resonance frequency should take 
place across the transition point between the two degen- 
erate states. This is at odds with an abrupt onset of 
energy absorption at B > £>* that is expected in the 
absence of the Kondo coupling. 



APPENDIX A: THE SCHRIEFFER- WOLFF 
TRANFORMATION. 

In this Appendix we report the details of the Schrieffer- 
Wolff transformation, leading to the Effective Kondo-like 
Hamiltonian (fTcf). 

The Hamiltonian matrix elements Hij between states 
= 0: 1) 2, 3) defined in and below eq.@ are in- 
troduced in and below eq. ^ . The Hamiltonian operator 
(^) projected onto the subspace S generated by the two 
degenerate levels (a, (3 = 1, 2) is: 



H eff = (Hi + e d )+J2M af} 

a/3 

The operators M a p are given by: 



(Al) 
(A2) 



The basic approximations we introduce in order to 
compute the M a p consists in keeping only lead excita- 
tions with energy around the Fermi level, whose energy 
is negligible. Moreover, we approximate E at the denom- 
inators with e d . Thus, the matrix elements become: 



M a t 



£d U + e d 

Their explicit form is: 

Moo = —d) a d a '^2v* aQ {q)v aQ (q )(() q cr a 4>\, 



(A3) 
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Maa = d^da } 



'I 'I 



v a3 (q) v &3(q') v a o(q')v a0 (q) 



U + e d 



e<2 



q&a r 1 cr Q 



(A4) 



Here is a ^ a and <Ji(2) =1 (?)■ 

The Effective Hamiltonian can be expressed as the fol- 
lowing operator: 

H e " =H l +e d + Y} A{ -1) + B(q)S z } + 



</</ 



'/'/ 
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J±(«,g , )[5-^ T ^ i + 0j^ g / T fi^] 



5* = S'qff is defined in the body of the paper, after 
eq. (|l^) . An enormous exemplification happens in eq. ( |A5| ) 
if one takes the amplitudes v to be independent on q and 
real. If this is the case, the coefficients are: 



A = l <± ■ 
2 e d 



20 . 



"23 



B = 



J 13 



J 20 
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e<2 



,,2 1 



f/ + £d 



Nevertheless, even in the non symmetrical case the po- 
tential scattering terms do not matter for the Kondo 
physics anyway, because they can be re-absorbed in a 
redefinition of the energy levels of the leads. Moreover, a 
"poor-man" 's scaling argument shows that the coupling 
J docs not scale as one lowers the energy cutoff (as 
confirmed by the corresponding differential equation in 
p7[ , as well). Hence, one can safely neglect it, as we 
(A5) have done in the paper. Finally, in the range of the pa- 
rameters relevant for our system, the fixed point of the 
unisotropic Kondo Hamiltonian we have obtained is the 
same as the fixed point of its isotropic limit, so that our 
model Hamiltonian ([n]) describes all the relevant physics 
of the problem. 

The results we obtained with the Schrieffer- Wolff 
transformation basically agree with the ones in [ff7| . 
(A6) However, at odds with jL7|, in working out the trans- 
formation we choose to keep the couplings to the N = 1 
and the N = 3 levels different. Indeed, both should be 
taken into account, unless one makes a special choice of 
V g , which appears to be unjustified [{33[. 



J 13 



"23 



U + e d 



APPENDIX B: RELEVANT EXPECTATION 
VALUES. 



J z = 



U + e d 



J± = 2 



^13^23 V W V 2 



U + e d 



J = 



J 13 



23 



U + e d e d 
and the Hamiltonian finally becomes: 

H eff =H l +e d + A + BS Z + U c p{0) + U s a(0) 



In this appendix we review the basic rules for the cal- 
culations leading to the results of Section IV and Section 
V. 

We start by computing the norm of the state \g x ) de- 
fined in eq.(p^),with the projector P g given by eq.(|T9|). 



N[g] = (g x \g x ) - (FS 7 X\(P g ) 2 \FS,X) 

= 2(1 + e) 



(Bl) 



where £ = \{\ — fff), as always. In order to reduce 
higher powers of the angular momentum in (P g ) 2 we 
make use of the general property of the spin-1/2 matrices, 
S a S b = + %e abc S c and of the following identities: 



+ {J z a z (0)S z + J x (a x (0)S x + a y {0)S y )J' p(0)S z } 

(A7) 

where p(0) = 4>\ a 4> q > a and 

The extra terms in the Hamiltonian are: i) a renor- 
malization of the relative position of the two degenerate 
levels, (this can be eliminated by re-tuning B) ; ii) a po- 
tential scattering term in the charge channel; Hi) a poten- 
tial scattering term in the spin channel; iv) an anisotropic 
Kondo coupling; and, finally, v) a spin-charge coupling. 

If the couplings to the states 1 and 2 are symmetrical, 
v a o and i>3 Q do not depend on a, the terms Hi) and v) are 
zero and the model becomes isotropic ( J z = J±_). This 
leads to eq.(p^) and to the matrix elements of eq. (|TT|) . 



aa z (0) = -(an ] (0)-an l (0)) 
(a*(0)) 2 = j(an T (0) + an^O) - 2 a 2 n T (0)n x (0)) 

H0)) 4 ^(a?(0)) 2 (B2) 



Using eq.(B2 ) we get the final result 



1 fi 

{P g f = (2£) 2 + - (1 - e) (^(0)) 2 - 8aa(0) ■ S (B3) 

Taking the expectation value of 5* on |A) and because 
(FS\an a (0)\FS) = \, (where an a (0) = j? J2 kk , c\ a c k , a ) 
we use 



(FS\aa z (0)\FS) =0; (FS\(aa(0)f\FS) 



(B4) 
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to obtain eq.(Bl). 

Next, we calculate the expectation value of the Kondo 
interaction Hamiltonian Hk denned in eg. (|T^) . Because 
Hk and P g commute, we get: 



E K ■ N[$ = J(FS, X\P 2 aa(0) ■ S\FS, A) 



2 J (B5) 



The interaction between the delocalized conduction elec- 
trons and the impurity provides corrections to the ground 
state energy of the FS, (FS\Hi\FS) = ~E GS = -N t D, 
that is again of order 0(1) in the number of particles. 
Such a correction is expressed as: 

AE C ■ N[£] = (FS, X\P g [Hi,P g ]\FS, A) = 

Ul-2O(FS,X\P g [H h (aa(0)) 2 }\FS,X) 



-4(FS,X\P g [H h aa{0)-S]\FS, A) 



(B6) 
Because 



We now calculate the r.h.s of eq.(| 
[Hi,an a (0)] = ^J2 qq ' £q{ c \ a Cq> a ~ c l' a c qa ), 
{FS, X\[Hi,an a (0)]\FS, X) = and, using Wick theorem 



we find 



(FS,A|ona(0)[Fi,ana(0)]|.F£,A) = 



Eqs 
2Ni 



(B7) 



Next, one more contribution to the kinetic energy is: 
(FS,X\aa a (0)[Hi,aa b (0)}\FS,X) = 



E, 



GS 



Ti- 



er" a" 



E, 



GS 



2Ni 



Finally, the variational estimate of the energy correction 
due to the coupling to the impurity, which is 0(1) in the 
particle number, is: 



1 



e(£,J) = -(AE c + E K ) 



2(1 -tj+e)-p/p 



(B8) 



Because J/D = u(0)J/Ni, we obtain eq.(pp[) of the text. 

In Section V we extend eq.(BJ ) to the SB ^ 0-case. 
The new quantity that appears in the problem is h\, 
defined by eq. (|24|) . All the operator identities we have 
proved so far still hold in the SB ^ 0-case. However, 
being h\ ^ 0, the average values of (ac?(0)) 2 is now given 
by: 

(FS,SB\(aa(0)) 2 \FS,SB) = jj [l + (h x ) 2 ] (B9) 

which we used in order to obtain e[£, h\] in eq. (p5|) . 

Finally, the average value of the z-component of the 
S*jj at nonzero h, which is required in eq.((2^) is given 
by: 



effl 



N[{, h x 



■(FS,X\P g S z P g \FS, A) 5B #o- 



(BIO) 



All other calculations in the text are straightforward. 
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FIG. 1. (a) Vertical Quantum Dot with contacts in a pillar 
configuration, (b) Grey-scale plot of the differential conduc- 
tance vs gate voltage V g and biasing voltage V s d (schematic), 
showing the CB regions (white) at various N. Operation point 
is marked by <g>. 
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FIG. 2. (a) General scheme of the energy level of the QD vs 
magnetic field B. Only few states are plotted (labels are de- 
fined in Section II). The level spacing of the parabolic lateral 
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action U = 3meV. Circled numbers represent the four dot 
states involved in the Kondo resonance, (b) Energy gain of 
the Kondo correlated state (not in scale) and removal of the 
degeneracy at SB — B ~ B* ( h oc SB /Tk'- see Section 
V). 
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Kondo conduction has been observed in a quantum dot with an even number of electrons at the 
Triplet- Singlet degeneracy point produced by applying a small magnetic field B orthogonal to the dot 
plane. At a much larger field B ~ B* , orbital effects induce the reversed transition from the Singlet 
to the Triplet state. We study the newly proposed Kondo behavior at this point. Here the Zeeman 
spin splitting cannot be neglected, what changes the nature of the Kondo coupling. On grounds of 
exact diagonalization results in a dot with cylindrical symmetry, we show that, at odds with what 
happens at the other crossing point, close to B„, orbital and spin degrees of freedom are "locked 
together", so that the Kondo coupling involves a fictitious spin 1/2 only, which is fully compensated 
by conduction electrons under suitable conditions. In this sense, spin at the dot is fractionalized. We 
derive the scaling equation of the system by means of a nonperturbative variational approach. The 
approach is extended to the B ^ B*-case and the residual magnetization on the dot is discussed. 

PACS numbers: 71.10.Ay, 72.15.Qm, 73.23.-b, 73.23.Hk, 79.60.Jv, 73.61.-r 



I. INTRODUCTION 

Low temperature transport properties of quantum dots 
(QD) are primarily determined by electron-electron in- 
teractions [1]. When a QD is weakly coupled to the con- 
tacts its charging energy is larger than the thermal en- 
ergy. Therefore, in the absence of a source-drain voltage 
bias, V s d, electrons can be added to the dot in a con- 
trolled way by changing an applied gate voltage, V g [2,3]. 
Conductance peaks are separated by Coulomb Blockade 
(CB) regions, in which the electron number TV is fixed 
(Coulomb oscillations). 

However, there is a well established experimental ev- 
idence for zero-bias anomaly in the conductance in a 
CB region when TV is odd [4,5]. This occurs when the 
temperature is further decreased and the strength of the 
coupling to the contacts increased. Such a behavior is 
attributed to the formation of a Kondo-like hybridized 
state between the QD and the contacts [6,7]. Kondo 
effect is not expected when TV is even, because of lack 
of spin degeneracy as the QD is in a singlet state. An 
exception is provided when a QD with TV even is in a 
triplet state because of Hund's rule. In this case the 
QD shows a very peculiar spin-1 Kondo Effect [8-11]. 
The Kondo anomaly is enhanced if the singlet state is 
close to be degenerate with the triplet. Such a regime 
can be achieved experimentally by applying a weak mag- 
netic field orthogonal to the dot, B±, because Hund's 
rule breaks down quite soon. Correspondingly, the QD 
undergoes a triplet-singlet transition, which eventually 
suppresses the Kondo effect. In the following we shall 
refer to the transition above as "TS crossing" . 

If the magnetic field is further increased, orbital ef- 



fects induce the reverse transition from the singlet to the 
triplet state at some "critical field" B± = 73* , in the gen- 
eral landscape of transitions to higher spin states [12,13]. 
In this paper we focus on this transition, in the following 
referred to as "ST crossing" . The accidental degeneracy 
between states with different total spin can restore the 
Kondo effect. However, as we will show below, the pres- 
ence of a substantial Zeeman splin splitting (Zss) makes 
the ST crossing very different from the TS crossing. Our 
analysis is an extension of the results of [14]. 

We consider a vertical QD with contacts located at 
the top and at the bottom of a pillar structure [2] (see 
Fig. la). Because the confining potential is chosen to be 
parabolic in the radial direction, single particle states on 
the dot are labeled by (n,m,a) where n is an integer, 
to is the orbital angular momentum and a the electron 
spin. We have chosen TV = 2 for sake of simplicity, but 
we believe that the pattern we describe is quite frequent. 
For example, it can be extended to the case when TV=4. 

Because of the Zss, the crossing occurring at the ST 
point involves the singlet and the S z = 1-component of 
the triplet state only. The analysis of the quantum num- 
bers of the TV ± 1 states visited by virtual transitions (co- 
tunneling processes ) shows that orbital and spin quan- 
tum numbers of the electrons involved are "locked to- 
gether" (Fig. 2a). Tunneling from the contacts does not 
conserve n, but, since the geometry is cylindrical, it con- 
serves to and the spin. In particular, if the dot is in the 
triplet state, only (to = 0, J,)-electrons can enter it, while 
is it is in the singlet state, then only (m = 1, t)-electrons 
enter. This envisages a one-channel spin 1/2-like Kondo 
coupling different from what occurs at the TS point. 

The existence of a Kondo effect between states belong- 
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ing to different representations of the total spin is very 
peculiar. Being the total dot spin of both the degener- 
ate levels integer (S = 0,1), spin compensation has to 
be incomplete. In this work we show that it is as if the 
total spin were decomposed into two fictitious spins \ 
which we refer to as S r and S e g [15]. We stress that 
these dynamical variables involve both orbital and spin 
degrees of freedom which are locked together. The mag- 
netic field _B* favors antiferromagnetic (AF) coupling be- 
tween the delocalized electrons of the contacts and S e g, 
which quenches it in the strong coupling limit. S r is the 
residual spin \ at the dot [16]. CB pins TV to an even 
value, but the correlated state that sets in will have half- 
odd spin, i.e., "fractionalization of the spin" at the QD 
site will take place [14]. B* only affects the magnetiza- 
tion related to S r , but we will not discuss such a point in 
this paper. On the contrary, detuning B from B* corre- 
sponds to an effective magnetic field SB = B — B*, which 
affects Kondo screening. 

On the other hand, at temperatures T > Tk ( Tk 
being the Kondo temperature), there is no Kondo corre- 
lation. Furthermore, being the dot in CB, the screening 
by the delocalized electrons is quite negligible. As long 
as B < B*, the total spin is 0. For B > B*, the total spin 
is 1. The spin density on the dot has a drastic change 
at B*, while the charge density is only slightly affected. 
This is shown in Fig. 3a: for B < B* (singlet state) the 
spin density is zero, while it is half the charge density for 
B > B* (triplet state S = 1 with z— component S z = 1). 

This shows that the total spin of the dot is the relevant 
dynamical variable in the correlation process that sets in 
when temperature is lowered, while the charge degree of 
freedom is frozen for B ~ 

In Section IV we set up a variational technique similar 
in spirit to Yosida's approach [18] for the construction of 
a trial wavefunction of the correlated Ground State (GS). 
We start from a Fermi sea (FS) of delocalized electrons of 
the contacts, times the impurity spin wavefunction and 
we project out states which do not form a singlet at the 
dot site. This is a modified version of the Gutzwiller 
projector and the results we obtain qualitatively repro- 
duce the strongly-coupled regime of the Kondo model 
[7,19-21]. Our formalism also allows us to discuss the 
SB ^ 0-case and to calculate "off-critical" quantities, like 
the magnetization, as a function of SB. In [14] we esti- 
mated the Kondo temperature Tk- It is slightly below 
the border of what is reachable nowadays within a trans- 
port measurement (Tk ~ tens of mK). Instead, an Elec- 
tron Paramagnetic Resonance (EPR) experiment could 
provide evidence for the fractionalization of the spin. At 
T < Tk, the system goes across the Kondo state when 
B is moved across _B* . Partial compensation of the dot 
spin due to the local screening implied by the Abrikosov- 
Suhl resonance takes place at the QD site. The residual 
magnetization is the one of a spin-1/2 in all the range 
B ~ B*. Therefore, in lowering the temperature, the 



response to the microwave field should also appear all 
the way from B < B* to B > £?* and a temperature- 
dependent Knight shift of the frequency should occur. 
Moreover, this feature should be detectable, due to the 
enhanced susceptibility of the Kondo state. 

Differently from the scenario described here, the TS 
crossing has a nonsymmetric behavior across the degen- 
eracy point. Spin 1 Kondo effect should take place for B 
prior to the degeneracy point, It is renforced at the TS 
transition, but no Kondo effect at all should occur when 
the GS is a singlet. 

The paper is organized as follows: 

• In Section II we show that B % makes only four levels 
involved in the Kondo coupling, in analogy to the 
single impurity spin ^ Anderson model (AM). We 
discuss the strongly coupled regime of the system 
by briefly recalling the scaling theory of the AM and 
specify what we mean as spin fractionalization. 

• In Section III we employ the Schrieffer- Wolff trans- 
formation to map our system onto an Effective 
Kondo model. We show that the correspond- 
ing interaction between S e g and the spin of the 
electrons from the Fermi sea is antiferromagnetic 
(AF), which generates the strongly coupled Kondo 
regime. 

• In Section IV we construct a trial GS using a mod- 
ified version of the Gutzwiller projection method. 
We calculate the energy, which depends on the ex- 
change interaction: j. We show that j scales ac- 
cording to Anderson's poor man's scaling [22]. 

• In Section V we extend our technique to the SB ^ 
0-case. Because of the presence of SB, a second 
scaling parameter, h, arises. We compute the en- 
ergy and the magnetization as a function of h and 
discuss how this is related to the bare field SB. 

• In Section VI we report our final remarks and con- 
clusions. 



II. THE MODEL FOR THE DOT AT THE LEVEL 
CROSSING (B = B t ). 

Using exact diagonalization [13], two of us studied few 
interacting electrons in two dimensions, confined by a 
parabolic potential in the presence of a magnetic field B 
orthogonal to the dot disk (z— direction). We chose the 
energy separation between the levels of the confining po- 
tential fkod to be of the same size as the effective strength 
of the Coulomb repulsion U (~ 4meV), which makes the 
correlations among the electrons relevant. The magnetic 
field favors the increase of the total angular momentum 
in the z— direction, M, as well as of the total spin of the 
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dot, 5, in order to minimize the Coulomb repulsion. The 
^-dependence of the spectrum of the QD energy levels 
can be monitored via a linear transport measurement, 
by attaching electrodes to it in a pillar structure (Fig. 1) 
[4,5,9]. At special values B = B*, one can identify sev- 
eral level crossings between states with different M and 
5. 

In this paper we shall focus onto a QD with TV = 2 
electrons, but the features we discuss here are appropri- 
ate for dots with any even TV. The parameters of the 
device can be chosen such that T3* is rather large, which 
lifts the spin degeneracy, because of Zeeman spin split- 
ting [23] . We define B % as the field at which the GS of the 
isolated dot switches from singlet ( 2 5°; 5 = 0, M = 0) 
to the triplet state lowest in energy (the state S z = 1 of 
the three states 2 T Q \ ±l : S = 1, M = 1 and S z = 0, ±1). 

Only two TV = 2 states are primarily involved at the 
level crossing, namely, the 5 = singlet and the S = 
1,M = 1,S Z = 1 component of the triplet (see Fig. 2a). 
At the same B, the low lying TV = 1 and TV = 3-GS's 
are spin doublets. We shall denote them as 1 D° 1 and 
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3 D\ , respectively. Lifting the spin degeneracy makes the 

2 

S z = \ components lower in energy in both cases. If V g 
is tuned in the CB region between TV = 2 and TV = 3 
(Fig. lb), only the four states listed above are primarily 
relevant for the conduction. 

We now construct a model by taking the 1-particle 
state , 1 D\, as the "vacuum" of the truncated Hilbert 

2 

space H (which we denote by |0)). Let us take |0) to be 
the initial state of the QD. The dot can exchange one 
electron at a time with the contacts. A (q, m, t)-electron 
entering the QD generates a transition between |0) and 
2 T/. Alternatively, the lowest energy state available for a 
(q, m, J,)-electron entering the QD is the singlet 2 5$ (here 
q is the modulus of the momentum of the particles and 
to is the azimuthal quantum number, which is conserved 
by the interaction with the dot) . It is useful to define the 
operators d\ 2 which create the two TV = 2 many body 
states \a) (a = 1, 2), belonging to H, by acting on |0): 

|l) = | 2 5 °) = dl|0); |2) = | 2 T 1 1 )=4|0) (1) 

The operator d[/d\ adds a (m = 0,i)/(m = l,f) elec- 
tron to the dot. Finally, the state |3) is defined as: 
1 3) = | 3 D\) = d\d\ 1 0). The operators d 12 defined 

2 

above cannot be associated to any single particle func- 
tion on the dot. Nevertheless, as far as low-energy exci- 
tations only are involved, it is a good approximation to 
take them to be anticommuting. To show this we refer to 
Fig. 4, where we report the calculated quasiparticle spec- 
tral weight, when an extra electron is added to the dot 
in one of the two possible GS's at B = B* and TV = 2. In 
Fig. 4a we show the weight for addition of an extra par- 
ticle to the state | 2 T/) to give a state with total angular 
momentum M = 1. This is close to 1 for a L-spin but 



is practically zero at low energies for an j' spin because 
states with 5 = 3/2 are higher in energies. According 
to the definitions of eq.(l), this implies that, while the 
state djd||0) can be normalized, one has (d|) 2 |0) « 0. 
Analogously, Fig. 4b shows that there is an energy shift 
of the spectral weight peak when a j' or a J, spin particle 
is added to the state | 2 5q). This shift corresponds to 
the Zeeman spin splitting between the states 3 D\ and 

2 

3 D 1 _ L . We ignore the latter that is placed at higher en- 

2 

ergy. In Section III we argue that inclusion of this state 
does not qualitatively change our results. Then, within 
this approximation that includes the lowest energy state 
only, we take (d|) 2 |0) £s 0. So, we are alleged to assume 
(d\ 2 ) 2 = 0. Hence, the d-operators behave as they were 
single-particle fermion operators, although, in fact, they 
create many-particle states. 

We write the Hamiltonian for the QD as: 

H D = e d (rai + n 2 ) + Un x n 2 (2) 

where n a = d^ a d a , e d is the energy of the degenerate 
TV = 2 levels and U is the charging energy for adding a 
third electron. 

We assume the electrons in the leads to be noninter- 
acting. Let bj qa be the annihilation operators for one 
electron in the left (j = L) and the right (j = R) con- 
tact, respectively, and let e q be the corresponding energy 
(independent of j and a for simplicity) . The Hamiltonian 
for the leads is: 

^ = EE zA° h ^ (3) 

qtj j=L,R 

We denote by Ho = He, + Hi the sum of the two Hamil- 
tonian terms given by eqs.(2,3). They do not change the 
particle number TV on the dot. Instead, TV is changed 
by tunneling between the leads and the dot, according to 
the above mentioned selection rules. Let T# and be 
the amplitudes for single-electron tunneling between the 
dot and the R or L contact, respectively. Removal of one 
electron from the contacts is described by the following 
linear combination of the 6j gcr 's: 

^ = At ^ , m d TLbL l° + T R h R<l°) ( 4 ) 
V I 1 il + \ l R\ 

For a point-like QD, there will be no tunneling term 
involving the combination of bjka operators linearly in- 
dependent of (f) qtT . The terms to be added are (a = 1,2): 

V a0 = X^«o(9)4<?W (5) 
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for TV = 1 -> TV = 2 and: 

V 3a = ^2 v 3a(q)dl(t> q a & , (6) 

Q 
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for TV = 2 — > N = 3(<5 ^ a) . Because spin is conserved 
in the tunneling is o\ =\., and o<i =\ and a ^ a. Her- 
mitian conjugate operators of V a o and V3 a decrease the 
electron number on the dot. 

This model Hamiltonian has the form of the nondegen- 
erate Anderson Model (AM). In [14] we have considered 
a special choice of the V g fine tuning (see Fig. lb) , which 
corresponds to the symmetric case: the singly occupied 
localized impurity level has energy = —U/2, while the 
doubly occupied one has energy + U = U/2. Also, 
Tj (j = L, R) have been taken independent on j. 

Integration over the contact fields <\> qa in the partition 
function leads to the constrained dynamics of the dot be- 
tween the two degenerate states \a), labeled by a = 1,2 
and coupled to an external fluctuating field X(t) [24]. 
As shown in [14], the main contributions to the partition 
function are the quantum fluctuations of an effective fic- 
titious spin i S e ff, due to the coupling among the dot 
and the leads. This spin flips from 1 to 2 because of 
cotunneling processes in which an (to = l,t) electron is 
added/removed in the dot while an (to = 0,1) electron 
is removed/added (see Fig. 2a ) Its quantum dynam- 
ics produces a correlated state between dot and leads in 
which the charge excitations, being higher in energy, are 
fully decoupled and suppressed. Four states are available 
prior to correlation, corresponding to S = 0,1, two of 
which are degenerate. They factorize according to the 
decomposition law: 



eff- 



(7) 



where S r is a residual spin 1/2 whose wavefunction fac- 
torizes close to B* . Correlation with the contacts is gov- 
erned by a single channel spin \ AF Kondo Hamiltonian 
for S e ff, which is derived in Section III and discussed in 
detail in Appendix A. In the corresponding GS, S e g is 
fully compensated as in the usual Kondo Effect at zero 
magnetic field. 

In Sect. Ill and IV we show that, provided that the 
coupling to the leads is symmetrical for a=l,2, the GS 
is doubly degenerate at B = £?* (the limiting states are 
derived from the \g x ) (A = ±) of eq.(16), having en- 
ergy e(£,j) given by eq.(20)). The degeneracy disap- 
pears when is SB = B — ^ (see Fig. 2b). Such 
a spin fractionalization closely resembles quantum num- 
ber fractionalization in strongly-correlated 1-dimensional 
electron systems [25]. 



III. EFFECTIVE AF KONDO INTERACTION. 

In the previous Section we have introduced the four 
states \i), i = 0,1,2,3 , which are primarily relevant to 
the low-temperature dynamics of the QD (see eq.(l)). 
Tuning V g as described in Section II (see Fig. lb) allows 
only the N = 2 states, |1) and |2), to be the initial and 



final states of any process. Let H be the subspace they 
span. The N ^ 2-states, |0) and |3), play a role only 
as virtual states allowing for higher-order tunneling pro- 
cesses. Then, it is desirable to describe the dynamics of 
the dot in terms of an effective Hamiltonian acting on H 
only. This is accomplished in this Section, by employing 
the Schrieffer- Wolff (SW) transformation at second order 
in the tunneling amplitudes. 

According to eqs.(2,3 ) the elements of "Ho are diagonal 
in the Fock space and are given by: 

H 00 = Hi; Hn = H22 = Hi + Ed\ i?33 = Hi + U + 2ej, 

(8) 

#03 = H12 = 0, while the tunneling part of the Hamil- 
tonian of eq.s (5,6) provides the remaining off-diagonal 
terms Hij = 0,1,2,3). 

The effective Schrodinger equation projected onto H is: 

p[{n - e) - yt(i - p)(n - EyW^) = o (9) 

where P is the corresponding projection operator. 

For the sake of simplicity we shall take Vs a (q) and 
v a o(q) to be independent of a. Then, to second order 
in v, eq.(9 ) becomes (see Appendix A for details): 

qq 1 a- 

E^'E^Wy^'Y (10) 

qq' a 

where er^ A , are the spin \ matrices (a = x,y,z) and 
S e ff is a spin-i representation of the rotation within S, 
which is well defined, provided the constraint ni+n^ = 1 
is satisfied. Indeed, in this case, (S e ff) 2 = §(rii + ri2 — 
2n\ri2) = §. More explicitly: S z e g = \(d\d,2 — d[d{) and 

the ladder operator is S~g = d\ d<i . 
The interaction matrix elements are 



Qq 



+ 



vl(q)v3(q ) , v (q)v^(q 
e d + U -e q 



V 



€d + U - e q 



(11) 



corresponding to the potential scattering term and the 
Kondo coupling term, respectively. 

The potential scattering term Q qq < provides a one-body 
interaction in the "charge" channel. It is irrelevant for 
the physics of the Kondo effect, which is related to inter- 
actions in the "spin" channel [20], and we shall disregard 
it in the following. 

By taking q at the Fermi surface, the Kondo term in 
eq.(10) becomes: 



H K = H eS -H l ~ JSe ff -J2^xx' 



(12) 



gg 
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Tuning V g such that the chemical potential of the leads 
is the zero of the single-particle energies, the Kondo cou- 
pling strength in eq.(ll) takes the usual form: 



J = \V\ 



1 



1 



€,{ + U € d 



(13) 



€d ?s —U/2 implies J > 0, i.e., the interaction between 
the spin density of the Fermi sea at the location of the 
dot and the effective spin S e ff, described by eq.(12), is 
antiferromagnetic. Thus, the low-energy physics of the 
system will be controlled by the strongly coupled AF- 
S = i-Kondo fixed point, where a spin singlet takes place 
at the impurity site. Because of the way S is related to 
S e & (eq.(7 )), total compensation of S e g leads to a spin 
S = \ at the QD. 

In the next Section we will make use of such an equiv- 
alence to construct a formalism for the strongly coupled 
fixed point of our system. 



U, where higher-energy states are "projected out" by the 
interaction and the GP method works successfully [27]. 
Our approach is similar in spirit to Yosida's variational 
technique [18]. However, here we are mostly interested 
in the "macroscopic" variable <?(0). Hence, it is the only 
lead operator involved in the construction of the vari- 
ational state. This makes our technique much simpler 
than the one in [18], because our approach requires just 
one variational parameter for the trial state. 

Here we list the basic assumptions concerning the 
model Hamiltonian and the trial variational state. We in- 
troduce our approach for the simple spin-^ Kondo Hamil- 
tonian (which corresponds to B = B„). In Section V we 
generalize it to the SB ^ 0-case. To simplify the nota- 
tion, we drop the suffix e g from S e g all throughout this 
Section and the next one. 

a. The Hamiltonian is (see eq.(12)): 

H = Hi+ H K = Hi + Jaa(0) ■ S (14) 



IV. GUTZWILLER PROJECTION. 

In Section II we have shown that cotunneling at the 
ST point strongly involves just two among the four avail- 
able states of the isolated dot, which can be described 
in terms of S e ff, including both orbital and spin degrees 
of freedom locked together. Instead, the other relevant 
degrees of freedom of the dot can be lumped into another 
pseudospin variable S r (S r = \) which decouples below 
Tk, because the dot is tuned at CB (eq.(7)). According 
to Section III the dynamics of correlations involves the 
true spin density of the delocalized electrons at the dot 
site 8 (0). Therefore, the study of the magnetization in 
the correlated state requires the knowledge of <r(0) . 

In Section III we discuss the assumptions under which 
the isotropic single channel spin | AF Kondo model of 
eq.(12) describes the physics at the fixed point. The cor- 
related state is a Nozieres local Fermi Liquid (FL) [26], 
with a spin singlet at the origin. In our case there is a 
substantial magnetic field B ~ B* and a small detuning 
SB ^ is likely to occur. 

In this Section, we construct an approximated GS of 
the strongly coupled system, described by H e % (eq.(12)). 
In order to describe the scaling of S e g toward total com- 
pensation, we use a variational method based on the 
Gutzwiller projection (GP) technique. This qualitatively 
reproduces the main features of the correlated singlet 
state. Detuning B off B % allows for probing the QD 
magnetization and for studying the magnetic response 
in the correlated state. So, in Section V we extend our 
approach to the B ^ £?*-case. 

As J increases, states other than a singlet at the im- 
purity become higher in energy and can be "projected 
out" from the physical Hilbert space. This is quite sim- 
ilar to what happens in the 1-d Hubbard model at large 



The impurity is located at x = 0. The scattering 
is s-wave, so that the model is effectively 1-d and 
can be defined on a lattice with spacing a and pe- 
riodic boundary conditions at L = TV/a. (TV/ is the 
number of lattice sites ) [28] . Hk provides an effec- 
tive interaction among the delocalized electrons. At 
each stage of the scaling process, we take the Slater 
determinant \FS) = Y\ a Yl q <o c l<r\fy as a reference 
state (FS: Fermi sea). The c qa operators annihi- 
late the quasiparticles and |0) is the quasiparticle 
vacuum. The spectrum is linearized around the 
Fermi point qp. D is the bandwidth. Therefore, 
the lead Hamiltonian is Hi = ^2 qa . e qa c qcr c qcr , where 
e qa = VFq- Here vf is the Fermi velocity and q is 
the momentum measured with respect to the Fermi 
momentum. The bandwidth is 2D = 2ttvf/<i- The 
density of states at the Fermi level is v(0) = TV/ /D 
and is assumed to be constant during the scaling. 

b. One can represent the real space field operator in 
terms of the c-operators: 

It satisfies the anticommutation relations 
m^^)}= 1 J a yS x y. 

The quantity an a (0) = aipl.(0)ip a (0) plays the 
role of the occupation number at the impurity 
site. This is consistent with the operator identity 
(an a (0)) 2 = an a (0). The corresponding spin am- 
plitude due to the delocalized electrons at x = is 
a<r(0) = ^a - 0a(O)<? aa ''0a' (0)- The average value of 
(aa(0)) 2 on \FS) is: (FS\(aa(0)) 2 \FS)=3/8. This 
is a consequence of averaging over configurations 
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with zero or double occupancy at x = (total spin 
0) and configurations with single occupancy at the 
same point (total spin 1/2). 

The correlation between the impurity and the FL is 
accounted for by projecting out of the uncorrelated 
state the components with a triplet or a doublet of 
the total spin at the impurity site, S + aer(0). 

Denoting the projector as P g , the variational trial 
state is defined as: 

\g x )=P g \FS,\)=P g (\FS)<E>\\)) . (16) 

Here | A) is an eigenstate of S z with eigenvalue ^A = 
±i and g is the variational parameter. The explicit 
form of the Gutzwiller projection operator P g will 
be given below. So far, it is enough to say that all 
the components of the state other than a localized 
singlet at x = are projected out, as g changes 
from zero to the limiting value, g = 4/3. 

At B = B* , we assume the usual one parameter scal- 
ing process when D — > D — SD. g approaches 
4/3, and, in order to to guarantee that at each 
step the total energy is at a minimum, the pa- 
rameter v(Q) J is correspondingly renormalized. If 
SB = B — B* ^ 0, a second parameter appears in 
the model, related to SB. We shall refer to such a 
parameter as h and will introduce it in the next Sec- 
tion. It is generated by the shift in the quasiparticle 
energies, according to their spin and will determine 
the "off-critical" magnetization (see point e). 

According to eq.(7), the effective magnetization of the 
dot is given by the average of the total spin at x = 
on the trial state \g x ) (shortly denoted by {■■■)g 1 sB)'- 

lim {M d (g,SB)-M d (g,SB = 0)} 

3^4/3 

~ lim g*HB [{S z eff ) g ,SB + {aa z (0)) g , 6B ] 

g—y4/3 

= M imp + m(h) . (17) 

Here we have temporarily restored the suffix e ff, 
for clarity. The dot magnetization M4 includes the 
local moment on the impurity site Mj mp , and the 
induced magnetic moment of the delocalized elec- 
trons from the contacts, m(h). 

At SB = both M imp and m(h) vanish, due to the 
compensation of the effective spin S e g of the stan- 
dard spin i AF Kondo problem. Md(g,SB = 0) is 
the residual magnetization at B = £?* associated 
to the expectation value of S r of eq.(7) and is not 
discussed in this paper. However, the spin suscep- 
tibility of S r is expected to be much smaller than 
the Kondo susceptibility of S e g. 



One remark concerns the fact that the giromag- 
netic factor in eq.(17) might not be the bare one, 
appearing in the electron Zeeman spin splitting. 

In Section V we work out the form of M, mp and 
m(h) in the strongly coupled scaling regime. 

Let us now construct P g . 

According to Noziere's FL picture of the correlated 
Kondo state [26], the main effect of the interaction be- 
tween the FS and the impurity taken at the AF fixed 
point is to fully project out the components of the trial 
state other than a singlet of the total spin at the impurity 
site. Let us consider the operator 

aS(0) ■ S = ±[(a<?(0) + Sf - \ - (aa(0)) 2 } (18) 

It will give zero when acting on a state that has a doublet 
on the impurity site, (i.e. ((a<r(0)) 2 ) = 0), —3/4 when 
acting on a singlet, 1/4 when acting on a triplet. Among 
the three possibilities listed above, the state correspond- 
ing to the triplet is the one highest in energy. We fully 
exclude it from the trial state from the outset and grad- 
ually project out the doublet by varying the parameter g 
from to 4/3. Accordingly, we define P g as: 

P g =[\-\g- gaa(0) ■ (1 - 4aa(0) • S) = 

(l-^g + g(aa(0)) 2 ^-4aa(0)-S (19) 

This operator projects out the components other than 
the localized singlet. When g=4/3 the projection is com- 
plete. 

In this section we report the result for the zero mag- 
netic field case. The key points of the calculations are 
summarized in Appendix B. 

The expectation value of the Hamiltonian in eq.(14) 
differs from the energy of the reference Fermi sea \FS), 
E Q = —2TTVFNi/a, just by a term of order 0(1) in the 
particle number AT/. In units of D, this correction is: 

€{ ^ j) - d { m ~ Eo ) ~ 2 i+e — ( 0) 

where N[£] is norm of the trial state (see eq.(Bl)) and 
we have introduced the dimensionless variables j = 
3v(0)J/(8Ni) and £ = §(| — g) for convenience. 

It is worth stressing that e(£, j) does not depend on 
the polarization of the impurity, A. Such a degeneracy 
disappears when SB ^ [29]. 

For each j oc J/D, eq.(20 ) takes a minimum as a 
function of £ (i.e. g). The minimum £(j) w.r.t. j is given 
by: 

I(i) = -3 + VP + 1 (21) 
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Both the triplet and the doublet are fully projected out 
of the trial GS when g — > 4/3, i.e. when £ flows to zero. 
Eq.(21) shows that this limit corresponds to j — > oo. 

The derivative of eq.(21) w.r.t. £, reproduces the poor 
man's scaling law for j [22] provided g is small (that is, 
£ is close to 1/2) : 



dj 



d(lnO 



(22) 



Eq.(22) has been worked out in [22] within perturbation 
theory, provided g is small, that is, £ is close to \. In our 
simplified approach , eq.(21) (not eq.(22) !) holds all the 
way down to the fixed point. 



V. THE MAGNETIC MOMENT 

In this Section we qualitatively discuss what happens 
to the energy and to the magnetic moment of the dot 
when a magnetic field is present and B is possibly de- 
tuned off . 

Here we do not consider the in detail the magnetic mo- 
ment associated to S r . Instead, we focus on the one in 
the third line of eq.(17). We generalize the construction 
of Sect IV to the SB ^ 0-case by shifting the energies 
of the quasiparticles of the FS according to their spin. 
Two corrections to the magnetic moment arise due to 
this shift: a) (S e ff) gt B may not scale down to zero, b) an 
induced moment of the delocalized electrons m(ft) may 
arise. We first discuss the latter contribution. In the ab- 
sence of correlations, because B % is the degeneracy point, 
the energy difference of electrons involved in cotunneling 
at the Fermi level is the same for both spins. There- 
fore, provided also the density of states of the conduction 
electrons is the same at the Fermi energy for both spins, 
B % will only be responsible for changes in the tunnel- 
ing Hamiltonian due to orbital effects [30], which do not 
influence the flow to strong coupling in a relevant way. 

This is not the case if SB ^ 0. Qualitatively, there will 
be changes in the number of the conduction electrons 
involved, according to their spin a = ±i as follows: 



N a - N t /2 
Ni/2 



tSB 



Ni 



SB 

r 

D 



(23) 



The interaction modifies eq.(23) by scaling, so, we gen- 
eralize it by defining a dimensionless quantity ft in terms 
of the expectation value of the quasiparticle spin at the 
origin evaluated on the reference state \FS): 

(FS,\\aa z (0)\FS,\)sB^o = 
(FS, A|| (an t (0) - an 4 (0)) \FS, A) = \h x , (24) 

(h\ = A • ft, where A = ± is the polarization of the impu- 
rity ). Eq.(24) states that the magnetic field produces a 



spin density at the impurity site proportional to ft. When 
SB = 0, rotational symmetry of the FS implies ft = un- 
less a local paramagnetism takes place, so that ft is a 
homogeneous function of SB. 

Eq.(24) is enough to generalize eq.(20) to the SB ^ 0- 
case. Using the results of of Appendix B, we get: 



[2(i-£+a+Mi-20](i 



h x ) - 2j(l + h x ) 



[(i + h x )+ea-hx) 



(25) 



As in the previous Section, the minimum of eq.(25) 

w.r.t. £ {de/di = 0) fixes T(j,h) ~ (1 + h)/2j. The 
states |<? A )(A = ±) are no longer degenerate, except for 
ft = and they represent the ground and excited state, 
depending on the sign of ft ( Fig. 2b). At ft = the 
delocalized electrons do not have an intrinsic magnetic 
moment close to the fixed point: m(h = 0) = 0. Because 
of the degeneracy, the magnetization at ft ~ is obtained 
from the average E/J = (e + + e_)/2. Minimization with 
respect to £ yields the leading large j-correction to the 
energy for ft ^ 0, due to the delocalized electrons of the 
leads: 



J 4 V Jo Jo 



(26) 



Hence, the contribution from the delocalized electrons to 
the magnetic moment is linear in ft as expected from the 
exact result [20]. However, the relation between ft and SB 
cannot be established within this model, unless we state 
some direct link between the strong coupling fixed point 
and the uncorrelated state previous to scaling. This link 
can be inferred from the educated guess that the response 
of the conduction electrons to SB, given by eq.(24) in the 
absence of correlations, smoothly evolves into the one of 
the correlated state, m(ft), in the form: 

m(ft) = lim <<7 A |aa Z (0)lA> ~ \^fh x (27) 

(here Tk is the Kondo temperature). 

Indeed, by using the insight we have from the exact 
solution, we infer that, the fraction of delocalized elec- 
trons involved in the correlated state is expected to be 
0(Tk/J)- Because Tk is a scale invariant, both ft and j 
scale to larger values in the same way as £ — > 0. 

Then, by comparing eq.(27), supplemented with the 
derivative of eq.(26), to the magnetization of the uncor- 
related state, Tj(Nf - N±)/Ni = g* HbSB j '4D , we obtain: 



1 dE 
2J~dh 



SB_ 



(28) 



Thus, we recover the expected result that the spin sus- 
ceptibility is oc l/T K [20]. 

This concludes our discussion of m(h). 
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According to eq.(17), the total magnetization also in- 
cludes the term M, mp . This term just provides a sub- 
leading correction to the total magnetization. 

Indeed, M imp is given by 

1 £ 2 

Mi mp (Z,h x ) = -n B g*\ 1+h * - , (29) 

l-hx + ? 

i.e. this correction to the fractional spin of the dot de- 
creases as the second power of £ when £ decreases in the 
strong coupling limit. 

VI. CONCLUDING REMARKS 

Kondo coupling is the striking realization of non- 
perturbative tunneling between a QD in the CB region 
and contacts. It requires the degeneracy of the GS dot 
level, what certainly takes place when TV is odd. In the 
even- TV case, the GS can be a triplet, because of Hund's 
rule and corresponding partial screening of the 5 = 1 
spin on the dot has been observed [8,9]. The proximity 
to a Triplet-Singlet (T-S) transition enhances the cou- 
pling and increases the Kondo temperature [10]. A small 
magnetic field driving the system across the level crossing 
can tune the intensity of the effect due to the T-S con- 
version. This enhanced Kondo conductance was indeed 
measured for TV = 6. 

In [14,17] a different mechanism for the Kondo effect 
in an even-TV QD has been proposed, considering states 
belonging to different representations of the total spin, 
5=0,1, at the accidental crossing between the singlet and 
the triplet state lowest in energy. In [17] B is taken in 
the plane of the dot and the crossing is attributed to the 
Zeeman spin splitting term. Because the latter is usually 
very small, it is unlikely that this can occur anywhere else 
but at the point mentioned above. Then the mechanism 
discussed by [8,10] seems to be more likely [31]. 

On the contrary, by increasing B± up to a suitable 
value £?* , the degeneracy of the triplet state is lifted and 
the singlet state ceases to be the GS of the system. In [14] 
we estimate quantitatively the parameter values of the 
optimal device by using results of exact diagonalization 
for few interacting electrons in a dot [23] . A vertical QD 
with a magnetic field along the z— axis, undergoes a S-T 
transition mainly due to the Zeeman term, which favors 
larger angular momenta along z, and to the correlation. 
The spins of the electrons in the contacts are still unpo- 
larized at the Fermi energy if the contacts are bulk (i.e. 
3-dimensional) metals. This allows Kondo coupling to oc- 
cur at the crossing between the two lowest levels among 
the four with 5 = 0, 1. In [17] the zero voltage anomaly in 
the differential conductance is discussed by including just 
these two levels. The fact that two different spin states 
are involved is reflected in the asymmetry in the trans- 
mission probabilities w.r.t. the spin index. This can be 



traced back to a term which couples the charge and spin 
degrees of freedom in the exchange Hamiltonian that can 
be derived from a Schrieffer- Wolff transformation. How- 
ever, this term is irrelevant as the system flows toward 
the strongly-coupled fixed point, consistently with the 
full decoupling of the charge and the spin degrees of free- 
dom in such a limit, due to the widely different energy 
scales for the two excitations [20]. The Schrieffer- Wolff 
transformation (see Section III) maps the problem onto 
the AF spin \— Kondo Hamiltonian. Coupling is better 
described by splitting the total dot spin S = 0, 1 into two 
spins i, which we refer to as S r and S e g (see eq.(7)). S e g 
is totally screened, independently on , while a residual 
spin i survives, whose magnetic moment depends on 

In this work we have concentrated on the magneti- 
zation at the strong coupling fixed point. A variational 
approach is used for the GS, based on the Gutzwiller pro- 
jection method, which, as far as we know, has never been 
applied to the Kondo problem. This technique allows for 
studying a small detuning of B off B % and for discussing 
the dependence of the magnetic moment on SB = B — B*. 
The basic assumption is that SB ^ does not move the 
system away from the strongly coupled state. The scale 
of magnetic field at which the Kondo state is disrupted, is 
expected to be quite large (it was estimated in [14] to be 
of the order of 1 Tesla). A similar conjecture in the case 
of magnetic impurities in metals has been formulated by 
Y. Ovchinnikov et al. [32], starting from a perturbative 
approach in a reduced Hilbert space. 

When SB ^ 0, besides the Kondo coupling parameter 
j oc i/(0) J, our variational correlated state includes a 
new coupling parameter which is related to the detuning 
magnetic field, h\ oc SB/Tk- As the system flows toward 
strong coupling, this parameter increases together with 
j. Within our formalism we show that j scales according 
to Anderson's poor man's scaling, provided the ration 
h/j keeps small (< 1). 

We have proved that the magnetic moment on the dot 
has a term linear in h (see eq.s(17,27)), which arises 
from the screening due to the electron density of the 
Abrikosov-Suhl resonance located at the dot site. Still, 
more work has to be done in order to provide the relation 
between h and SB (in our formalism h is introduced as 
the shift of the quasiparticle energies corresponding to 
u(0)SB). 

The QD GS includes what we referred to as a spinon, 
carrying a magnetic moment S r = 1/2 but no charge. 
The splitting in energy with magnetic field arising from 
this term will depend on the spinon wavefunction itself. 
At this stage, an estimate of the renormalized giromag- 
netic factor is impossible. However, the residual half in- 
teger spin on the dot, together with the screening effects 
in the GS at non zero SB, has important experimental 
implications. 

In [14] Tk has been estimated to be rather low, i.e., 
of the order of tens of mK. At such a low temperature a 
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transport measurement is cumbersome. We proposed to 
probe the fractional magnetization of the QD by means of 
a magnetic resonance experiment with microwaves. The 
energy splitting of the dot magnetic moment can be de- 
tected by observing resonant absorption as B is slightly 
detuned from B*. 

A Knight shift should be measured, depending on the 
sign of SB, which is a consequence of the term m(h) 
in eq.(17). In Fig. 3 we plot the radial behavior of the 
z— component of the total spin in the GS when B is 
slightly less or larger than B*. The dot is isolated and 
in the Coulomb blockade regime at TV = 2. In absence of 
Kondo coupling, screening from the delocalized electrons 
of the contacts does not occur and the magnetic mo- 
ment changes abruptly from zero to its maximum value 
at the transition. On the contrary, in the Kondo regime, 
the spin susceptibility of the conduction electrons on the 
dot is expected to be large and a temperature dependent 
screening should be present in the neighborhood of B % . 
This gives rise to a smooth crossover in the paramagnetic 
resonance close to T3* at low temperature. 

In conclusion, we have discussed the new possibility 
that Kondo effect arises in a QD with TV = 2 at low 
temperature, when the dot is tuned to the ST point and 
the coupling to the leads increases. This happens when 
the GS of the unperturbed dot changes from 2 Sq (B < 
-B*) to 2 Tl (B > due to level crossing. 

The theory of [14] and of this work is not applicable 
to the TS crossing. Orbital effects are dominant in our 
case in producing the crossing as well as the properties 
of the many body states with TV + 1 or TV — 1 particles, 
as described in the text. All our results about the nature 
of the states involved are totally independent on the Zss. 
However, at the ST point the Zeeman spin splitting is 
quite substantial, and favors a certain component of S z . 
Its role becomes important in two respects: 

a) Only two TV-electron levels cross at the TS point 
(not four oh them, as is the case of the TS point), which 
we label by 1 and 2. 

b) states at TV, TV + 1 and TV — 1, although they 
are many-body levels, have quantum numbers such that 
there is just one channel by which an electron can be 
virtually added or subtracted to the dot in the TV-state, 
because spin and orbital momentum of the extra electron 
are strictly "locked together" in the selection rules. 

The spin of the QD becomes \ without changing the 
average occupancy [14]. We propose to probe spin frac- 
tionalization with a magnetic resonance experiment. A 
continuous shift in the resonance frequency should take 
place across the transition point between the two degen- 
erate states. This is at odds with an abrupt onset of 
energy absorption at B > B % that is expected in the 
absence of the Kondo coupling. 



APPENDIX A: THE SCHRIEFFER- WOLFF 
TRANFORMATION. 

In this Appendix we report the details of the Schrieffer- 
Wolff transformation, leading to the Effective Kondo-like 
Hamiltonian (10). 

The Hamiltonian matrix elements Tffy between states 
N)i = 0j 1j 2, 3) defined in and below eq.(2) are in- 

troduced in and below eq.(8). The Hamiltonian operator 
(9) projected onto the subspace H generated by the two 
degenerate levels (a, (3 = 1,2) is: 



H e ff = (H l +e d )+Y / M a(3 
The operators M a p are given by: 
M »V = H °°Ehr i Hof> + H ^E-U-2e d -H l 



(Al) 

H 3 p 
(A2) 



The basic approximations we introduce in order to 
compute the Ai a p consists in keeping only lead excita- 
tions with energy around the Fermi level, whose energy 
is negligible. Moreover, we approximate E at the denom- 
inators with e d . Thus, the matrix elements become: 



P 



HaoHop H a3 H 3 p 



£d U + £d 

Their explicit form is: 

M aa = — d a d a ^v* a0 (q)v a0 (q )4>qoA\ 

QQ 

1 



(A3) 



U + e d 



QQ 



Maa = d a d & ^ 



v* a3 (q)va3(q') Vao(q')v* ao (q) 



U + e d 



b ] 6 • 
(A4) 



Here is a ^ a and cr 1(2) =1 (t)- 

The Effective Hamiltonian can be expressed as the fol- 
lowing operator: 

H e ff =H l+ e d + J2l A (v) + B(q)S z } + 



^ c M')[^ t <Vt + ^l<Vi] + 



E 



QQ 



Us{q,q')\[<t>\^ q ^ ~ <t>' qi <t> q ' i\ + 



At 
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J±(g,g')[5-0! t 



'9' 4 



+ 



(A5) 



5 = S e ff is defined in the body of the paper, after 
eq.(10). An enormous exemplification happens in eq.(A5) 
if one takes the amplitudes v to be independent on q and 
real. If this is the case, the coefficients are: 



A=- 
2 



U, 



l«io + "20 . B 



J 2Q 



J W 



viz + v 2 



13 



^ + e d 



+ 



vl + v 2 



20 



Ed 



(A6) 



Nevertheless, even in the non symmetrical case the po- 
tential scattering terms do not matter for the Kondo 
physics anyway, because they can be re-absorbed in a 
redefinition of the energy levels of the leads. Moreover, a 
"poor-man" 's scaling argument shows that the coupling 
J does not scale as one lowers the energy cutoff (as 
confirmed by the corresponding differential equation in 
[17], as well). Hence, one can safely neglect it, as we 
have done in the paper. Finally, in the range of the pa- 
rameters relevant for our system, the fixed point of the 
unisotropic Kondo Hamiltonian we have obtained is the 
same as the fixed point of its isotropic limit, so that our 
model Hamiltonian (10) describes all the relevant physics 
of the problem. 

The results we obtained with the Schrieffer- Wolff 
transformation basically agree with the ones in [17]. 
However, at odds with [17], in working out the trans- 
formation we choose to keep the couplings to the TV = 1 
and the TV = 3 levels different. Indeed, both should be 
taken into account, unless one makes a special choice of 
V g , which appears to be unjustified [33]. 



J 13 



"23 



U + e d 



+ 



u 20 



'10 



APPENDIX B: RELEVANT EXPECTATION 
VALUES. 



J z 



'13 



+ V. 



23 



'10 



+ v. 



20 



U + e d 



J± = 2 



^13^23 ^10^20 



U + e d 



J 



J 13 



J 23 



J 10 



J 20 



U + e d e d 
and the Hamiltonian finally becomes: 

H e f f =H l + e d + A + BS Z + U c p(0) + U s a(0) 



In this appendix we review the basic rules for the cal- 
culations leading to the results of Section IV and Section 
V. 

We start by computing the norm of the state \g x ) de- 
fined in eq. (16), with the projector P g given by eq.(19). 



N[g] = {g*\g*) = (FS,\\(P g y\FS,\) 

= 2(1+ a 



(Bl) 



where £ = \ — \g), as always. In order to reduce 
higher powers of the angular momentum in (P g ) 2 we 
make use of the general property of the spin-1/2 matrices, 
S a S b = s±_ + i e abc S c and of the f n owin g identities: 



+ {J z a z (0)S z + J ± (a x (0)S x + a y (0)S y )J' p(0)S z } 

(A7) 



E 



QQ 



qa^q a 



and 



where p(0) = 

<?(0) = 4>\ a °<*p4> q ' p 

The extra terms in the Hamiltonian are: i) a renor- 
malization of the relative position of the two degenerate 
levels, (this can be eliminated by re-tuning B) ; ii) a po- 
tential scattering term in the charge channel; Hi) a poten- 
tial scattering term in the spin channel; iv) an anisotropic 
Kondo coupling; and, finally, v) a spin-charge coupling. 

If the couplings to the states 1 and 2 are symmetrical, 
v a o and Vs a do not depend on a, the terms Hi) and v) are 
zero and the model becomes isotropic ( J z = J±). This 
leads to eq.(10) and to the matrix elements of eq.(ll). 



aa z (0) = -(an t (0) - an 4 (0)) 
(aa(0)) 2 = |(an t (0) +an ; (0) - 2a 2 n t (0)n ; (0)) 

(aa(0)) 4 = ^(a ( ?(0)) 2 (B2) 
Using eq.(B2 ) we get the final result 
(P g ) 2 = (20 2 + y (1 - £ 2 ) (a<r(0)) 2 - 8ad(0) ■ S (B3) 

Taking the expectation value of S on |A) and because 
{FS\an a (0)\FS) = \, (where an a (0) = jf J2 kk ' 4o c *'«) 
we use 



(FS\aa z (0)\FS) = 0; (FS\(aa(0)) 2 \FS) 



(B4) 
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to obtain eq.(Bl). 

Next, we calculate the expectation value of the Kondo 
interaction Hamiltonian Hk defined in eq.(12). Because 
Hk and P g commute, we get: 



E K ■ N[i] = J(FS, \\P*aa(0) ■ S\FS, A) 



"2 J ( B5 ) 



The interaction between the delocalized conduction elec- 
trons and the impurity provides corrections to the ground 
state energy of the FS, (FS\H t \FS) = -E GS = -NiD, 
that is again of order 0(1) in the number of particles. 
Such a correction is expressed as: 

AE C ■ N[£] = <FS,A|P s [#,,P fl ]|FS,A) = 

|(1 - 20(FS, \\P g [H t , (aa(0)) 2 ]\FS, A) 



-4(FS,X\P g [Hi,aS(0)-S\\FS,X) 



(B6) 



We now calculate the r.h.s of eq.(B6). Because 
[Hi,an a {0)] = jfT,gg'£q(4a c q ' a - c l' a c Q»), wefind 
{FS, \\[Hi,an a (0)]\FS, A) = and, using Wick theorem: 



(FS,A|an a (0)[ff ; ,an a (0)]|FS, A) 



Eqs 
27V 



(B7) 



Next, one more contribution to the kinetic energy is: 



{FS,\\aa a (0)[H l ,aa b {0)}\FS, A) 



E, 



GS 



Tr 



TT 



E, 



GS 



2iV ; 



Finally, the variational estimate of the energy correction 
due to the coupling to the impurity, which is 0(1) in the 
particle number, is: 



(AE C + E K ) 



2{i-i+e)-p/D 
(1+a 



(B8) 



Because J/D = v(0)J/Ni, we obtain eq.(20) of the text. 

In Section V we extend eq.(B8 ) to the SB ^ 0-case. 
The new quantity that appears in the problem is h\, 
defined by eq.(24). All the operator identities we have 
proved so far still hold in the SB ^ 0-case. However, 
being h\ ^ 0, the average values of (a<r(0)) 2 is now given 
by: 

{FS, SB\(aa(0)) 2 \FS, SB) = | [l + (/i A ) 2 ] (B9) 

which we used in order to obtain e[£, h\] in eq.(25). 

Finally, the average value of the ^-component of the 
S z e g at nonzero h, which is required in eq.(29) is given 
by: 



{Slff) 



{FS,\\P g S z P g \FS,\)sB?o = 

e{i-hi) 

4[{i + h x y + e{i-hi)] 



(BIO) 



All other calculations in the text are straightforward. 
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FIG. 1. (a) Vertical Quantum Dot with contacts in a pillar 
configuration, (b) Grey-scale plot of the differential conduc- 
tance vs gate voltage V g and biasing voltage V s d (schematic), 
showing the CB regions (white) at various N. Operation point 
is marked by ®. 
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FIG. 2. (a) General scheme of the energy level of the QD vs 
magnetic field B. Only few states are plotted (labels are de- 
fined in Section II). The level spacing of the parabolic lateral 
confinement of the QD is ftwo= 4meV , the Coulomb inter- 
action U = 3meV. Circled numbers represent the four dot 
states involved in the Kondo resonance, (b) Energy gain of 
the Kondo correlated state (not in scale) and removal of the 
degeneracy at SB = B — B*^0(/ioc SB/Tk- see Section 
V). 
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FIG. 3. Charge density as a function of the radius r. 
B < B„ (B > B„) should be understood as B slightly smaller 
(higher) than B*. The charge density p is only slightly af- 
fected when B goes across _B». The spin density a is zero 
when B < B„. It has the profile of the charge density as 
B > B, (a y = ip>). 
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FIG. 4. Spectral weights for addition of a spin up or down 
to the N=2 states, (a) Top panel-left: Spectral weight for 
addition of a (m = 1, f)-spin electron to 2 Tl; -right: Spectral 
weight for addition of a (m = 0, J.)-spin electron, (b) Bottom 
panel-left: Spectral weight for addition of a (m = l,f)-spin 
electron to 2 So; -right: Spectral weight for addition of a 
(m = 0, J.)-spin electron. Separation of the peaks is due to 
the Zeeman spin splitting 
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